Force on a moving point impurity due to quantum fluctuations in a Bose-Einstein condensate 
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An analytic expression is derived for a force on a weak point impurity arising from the scattering of quantum 
fluctuations in a slow-moving, weakly interacting, three-dimensional Bose-Einstein condensate at zero temper- 
ature. In an infinitely extended geometry, this force is shown to exist at any arbitrarily small flow velocity below 
Landau's critical velocity. Furthermore, this force is shown to be directly proportional to the flow speed. 

PACS numbers: 



(~~. ■ The exotic properties of superfiuid flow, especially the existence of a critical velocity below which the superfluid flows without 
, dissipation 1 1], continue to be the subject of active investigation after more than half a century. A recent letter |2] argued that 
the scattering of zero-temperature quantum fluctuations could contribute to a force that is present at all flow velocities, adding 
a new facet to this complex subject. Specifically, the effective critical velocity, defined by the minimum flow velocity at which 
O ■ an immersed object feels a force from the flow, was shown to be zero for a one-dimensional object in an infinitely extended 
superfluid modeled by a three-dimensional weakly-interacting Bose-Einstein condensate. However, the idealized 1-d potential 
used in this example has many peculiarities - such as the possibility of bound states, the presence of scattered waves at infinity, 
the force being negative and the changing sign of the effective mass - that potentially cloud the main idea that a drag force 
can be caused by quantum fluctuations at zero temperature. In this letter, we describe the calculation of a more experimentally 



realizable and theoretically transparent potential, namely that of a weak three-dimensional repulsive point impurity, which avoids 



the complications mentioned above. This calculation would be directly applicable to experiments involving trapped dilute BECs 
where untrapped atoms act as impurities We derive a simple analytical result that shows that the force on a stationary 
weak point impurity in the flow of an infinitely extended weakly interacting condensate is proportional to the flow speed and 
exists at any any arbitrarily small flow speed, and discuss how in a finite geometry this force may eventually cancel out due to 
£h ' backscattering effects. 

We consider a three-dimensional, weakly interacting condensate of density no, characterized by a delta-function interparticle 
contact pseudopotential whose coupling constant g is determined by the 2-particle positive condensate atom-atom scattering 
length a and the mass m of the atoms, i.e. g = 4irh 2 a/m, flowing at a speed c at infinity relative to a fixed weak impurity. 
Throughout this analysis, we assume the flow speed c to be less than the speed of sound. Although temperature is not well- 
defined in the scattering problem detailed in this paper, we use "T = 0" as a convenient description of the quantum state of the 
flowing condensate. While the scenario of a fixed impurity in a moving condensate is equivalent to that of a moving impurity 
in a fixed condensate, differing only by a Galilean transformation, we perform the calculation in the fixed frame of the impurity 
t— I ■ where the Hamiltonian is time independent. The potential of this fixed point impurity is described by r/S^ (r) where the coupling 
constant is given by n = 2nh 2 b/m, b being the scattering length that characterizes atom-impurity 2-body collisions; rj and b 
are defined as positive. We assume the condensate is dilute such that ^/noa :i <C 1 where no is the condensate density and the 
condensate's interaction with the impurity is weak. 

The general force arising from the flowing condensate on the fixed point impurity can be written in second quantized notation 

as 



; F = - / d 3 r(^( r )[V^( 3 )(r)]^(r)) T =o (D 

o 

where ip(r) and tp^(r) are field operators that describe the weakly interacting BEC flow and obey the standard boson commu- 
tation relations. The leading order contribution to the force is zero at flow speeds lower than the speed of sound; this result can 
be calculated assuming that the field operators are approximated by a macroscopic classical field representing the condensate 
whose behavior is governed in this mean-field approximation by the Gross-Pitaevskii equation (GPE) |4]. This absence of force 
below a critical velocity (in this case the speed of sound) verifies Landau's phenomenological picture, which does not explicitly 
include the quantum fluctuations. Based on this picture, Landau argued that the moving ground state would decay only above 
a critical speed (and only) by the emission of quasiparticles; otherwise the flow would be dissipationless [l]. Consequently, 
according to Landau's theory a stationary object in a superfluid flowing below this critical velocity would remain metastable and 
not experience any force. 

However, if we go beyond this mean-field picture by including quantum fluctuations, this conclusion that there is no force 
at low flow velocities no longer holds, at least when considering an infinitely extended medium. We briefly review the general 
framework for calculating the force (described in detail in |2|]). 

We approximate the Bose field operator that describes the flow by a large macroscopic classical field plus a small quantum 
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fluctuation operator (the small parameter in this case being given by the diluteness parameter ^noa 3 |5jfl), i.e. 

$(r) = *(r) + 4>(r) (2) 

where \&(r) is the macroscopic condensate that has been modified by the quantum fluctuations (see eq. [7). The force acting on 
an impurity can then be written as 

F = F fluc + F cond = r)V[{ft(r)4>(r)) T=0 + |*(r)| 2 ] r=0 (3) 

where the term F fi uc is the contribution to the force coming directly from the fluctuations and the term F con d is the contribution 
to the force coming from the condensate modified by the fluctuations. 

We expand <j) m terms of quasiparticle operators «k and a k , which by definition obey the standard boson commutation 
relations [dk, d k ,] = 5k, k' where <5 kjk ' is the Kronecker delta function and k is the dimensionless momentum normalized by the 
healing length. The quasiparticle operators are weighted by the quasiparticle amplitudes Uk(r) and Uk(r), i.e. 

<Xr) = £ («k(r)&k - <(r)4) , (4) 

k' 

where the sum is taken over all excited states and excludes the condensate mode. The normalization condition 
/ d 3 r (|uk(r)| 2 — |t>k(r)| 2 ) = 1 is imposed to force that the quasiparticle operators to obey the usual bosonic commutation 
relations. In order for the quasiparticle operators to diagonalize the weakly interacting Hamiltonian, the quasiparticle amplitudes 
must satisfy the two sets of differential equations known as the Bogoliubov-de Gennes equations |5l|6[: 

£u k (r) - * a «t(r) = E k u k (r) (5) 



£*« k (r) - (**) 2 M k (r) = -£!kVk(r), (6) 

where these equations are expressed in dimensionless variables with the length scale normalized by the healing length £ = 
{STrn^ay 1 / 2 ^ * denotes the complex conjugate, t = f + fj8^(r) - fj, + 2|*| 2 , f = -V 2 + y/2iq ■ V + q 2 /2, the 
dimensionless speed is given by q = c/c s , c s = y/n^g/m is the speed of sound, fi = 1 + q 2 /2 is the chemical potential 
determined by imposing f'°'(r) = 1 at r = oo, and fj — 2b/ an^ 3 . In this analysis, we assume that the point impurity 
is weak in the sense that fj <C 1. The energy eigenvalue associated with momentum state k for the moving BEC flow is 
Ek = \[2qk x + Eb where Eb = k\Jk 2 + 2 is the Bogoliubov dimensionless dispersion relation for a BEC at rest. 
The behavior of the condensate modified by the quantum fluctuations is given by the generalized GPE 1 8], 

(f + ^ 3 >(r) - M )*(r) + |*(r)| 2 *(r) + £ [2|« k (r)| 2 ^(r) - u fc (r)«£(r)¥(r)] - x(r)*(r) = 0. (7) 

k' 

The term proportional to 2^ k ' M k(r)wk(r) is ultraviolet divergent because of the contact potential approximation and must be 
renormalized (see for example |9]). The last term, %(r)^(r), ensures the orthogonality between the excited modes and the 
condensate |9], where %(r) = CkU k (r) where Ck = / d 3 r|^ , (r)| 2 [**(r)u k (r) + ^(r)wk(r)]- 

The effective scattering problem for the quasiparticle amplitudes determined by the Bogoliubov-de Gennes equations (eqns. 
|5J|6j 1 10] can be solved perturbatively with fj as the small parameter to give the total force on the impurity to leading order in fj: 

F x = <iV2Tr- 3 / 2 fj 2 pot 2 ^/noa 3 c J d 3 k(f fluc (k) + f cond (k)) (8) 

where ffi uc (k) and / C orid(k) are the k-dependent contributions of the force corresponding to F fi uc and F con d, respectively, and 
the zeroth order interaction pressure is given by po = gn\ /2. As we have assumed the flow to be only in the x-direction it follows 
that the direction of the force must be parallel to the x-axis. Henceforth, therefore, we drop the vector notation when referring to 
force. Note that because the leading order contribution to the force is proportional to fj 2 , it cannot be treated in a linear response 
formalism and a straightforward application of the fluctuation-dissipation theorem is not possible. The k-dependent contribution 
to the force arising directly from the quantum fluctuations is given by 



f fluc (k) = iV (k) / d 3 XX x [V]*(A, k) - Vi(A,k)] 



(9) 
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and the k-dependent contribution to the force from the modified condensate is given by 

iK. 



/ cond (k)-|d 3 A- 4 + 2A2 _ A ^ 2 
{x 2 [U (k) - 4V (k)] + V2\ x qU (k)} [V*(X, k) - Vi(A,k)] + V (k)[X 2 - V2X x q][U*(\,k) - #i(A,k)]) (10) 



where the tilde signifies the Fourier transform, i.e. f7i(r,k) = / d 3 Ae iA r f/i(A, k) and Vi(r,k) = / d 3 Ae iA r Vi(A ) k), and 
we have expanded to leading order in fj such that u^(r) = [Uq(\c) + f}U\(r, k)]e and v^(r) = [Vo(k) + fjVi(r, k)]e lk r . 



The zeroth order quantum amplitudes are given by Uo(k) = y \ \ ^Eb ^) aD ^ ^>(k) = y 5 y Eb — ^) w ^ tne ^ rsl order 

Fourier transform of the quantum amplitudes are given by J7i(A, k) and V\ (A, k). 

Below, we briefly sketch out the calculation for Ffi uc , noting that the same analysis can be applied to F con d, and then 
quote the results for both contributions to the force. Since Vi(A, k) is real, the only nonzero contribution to the force oc- 
curs when there is a pole on the real axis. Separating the integral over the Fourier variable A into Cartesian coordinates, 
we first perform the contour integration over X z . For simplicity, we separate out the characteristic equation of the cou- 
pled Bogoliubov equations, C(A,k) = -A 4 - 2A 2 (1 + k 2 + 2k x X x ) - AX x (k x + k 2 k x + kl\ x + E B q/V2 - X x q 2 /2), 
from Vi(A, k) = r„(A, k)/C(A, k) where r„(A, k) can be easily deduced from the Bogoliubov-de Gennes equations. 
The poles of r„(A, k) do not contribute to F con d because they are pure imaginary. Of the four X z poles of the charac- 



teristic equation, two of them, z%,2 = ±y — 1 — k 2 — 2k x X x — X 2 — X 2 — yl + Eg — 2EbX x + X x , are pure imaginary 
for all values of parameter space and thus do not contribute to F con d. The other two, Z3 4 = ±w(A^) 2 — X 2 where 



Xy = yj-l - k 2 - 2k x X x — X 2 + yl + Eg — 2EbX x + X 2 , can be real in a certain region of parameter space, namely where 
both X y is real and — X y < X y < X y . We perform the contour integral over X z to obtain 

ffluc(k) =4wV (k) [ dX x f3(X x ,k)L 3{X c )=0 (X x ) I V dX y 1 = (11) 



(Ag) 2 — X 2 

where /q(a<=)=o(A x ) is the indicator function which is defined as unity when Aa; is in subset S(A^) = and zero elsewhere; 

and (3(\ x , k) = X X Z3 ReSA s .=- X3 Vi(A, k), ReSj; =2 ;' being the residue at x = x' (note that f3{X x , k) is independent of X y as X y 
cancels out). The integral over X y gives a constant tt. Integrating over k, the integral over X x in the parameter regime where X y 
is real gives 



F fluc = -32V2ir 3 / 2 fj 2 Po eVnoa 3 sc / k 2 dkV (k) 







df k / + df k 0(X X ) 

-1 Jo Jf° J-xg 



(12) 



where k x = fc/fc, fj: = — qEg /\/2k(k 2 + 1); and, solving the equation A^ = for A^, we find A£ to be the only non-trivial 
real root. At large k, the term in brackets goes as k^ 1 and, since Vo(k) goes as k~~ 2 , the fc-integral has a logarithmic ultraviolet 
divergence and we can write out the leading order term as 

Ffi uc = ^j^TT 3/2 fj 2 P o^ 2 ^n Q a 3 sc qlnA uv (13) 

where A uv is the dimensionless ultraviolet cutoff 1 1 1 ]. Following a similar analysis, F con( i = 16-\/27r 3 / 2 7/ 2 po£ 2 \fnoa 3 c q In A uv . 
Thus, to leading order, the total force on an impurity in a moving, weakly interacting condensate at T = becomes 

F x - ( ^n 3 / 2 tf Po eV^o~aT c qlnAuv- (14) 

(Recall fj = 2b/ano^ 3 is assumed to be small, pa = gn^/2, £ = (87rnoa) _1//2 , and q is the Mach number, i.e. q = c/c s .) 
The ultraviolet divergence arises because of the singularity of the potential describing the point impurity. However, if we use 
a potential with a nonzero width in the x-direction, the result becomes finite. Therefore, in our problem, a sensible ultraviolet 
cutoff length is the dimensionless inverse of the scattering length describing the impurity-atom interaction 6, i.e. A uv w £/b. It 
is also reassuring to note that in the limit that the scattering lengths tend to zero, the force disappears. 

It is important to emphasize that this force directly proportional to the flow speed, reminiscent of Stokes' drag, exists at all 
flow velocities, including those below Landau's critical velocity. While we do expect this force to be dissipative - this can be 
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proven by computing the time dependence of the condensate energy in the lab frame, for example - the existence of this force is a 
necessary but not sufficient condition for the system to be dissipative. As discussed in 0O, this force is expected to eventually 
cancel out in a finite geometry due to backscattering effects. However, for the point impurity considered in the present letter, 
the time scale (a unique signature of this effect) over which the force cancels out should be on the order of the volume of the 
system divided by the product of the scattering cross-section of the impurity and the speed of sound. Moreover, as this force 
exists at all flow speeds in an infinitely extended geometry, one can speculate that the "persistence" of persistent currents arises 
from finite-size effects. 
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